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' Abstract. Dust plays crucial roles in galaxy formation and evolution. In 

V|j I the early epoch of galaxy evolution dust is only supplied by supernovae (SNe). 

With the aid of a new physical model of dust production by SNe, we constructed 
a model of dust emission from forming galaxies. We show the evolution of 
the spectral energy distribution (SED). Then we adopt this model to a local 
starbursting dwarf galaxy SBS 0335—052. Further wc discuss the SEDs of high 
redshift galaxies, and consider their observational feasibility. 
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en ' 1. Introduction 
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' Active star formation (SF) is followed by heavy element production from the 

\ birth and death of stars. Since the produced heavy elements generally exist in 

the form of dust grains, the dust grains absorb ultraviolet (UV) light and re-emit 
it in the far infrared (FIR). 

Indeed, there is an extreme category of high-z galaxies which have large 
amount of dust and are extremely luminous in the FIR and submillimeter 
^ fsubmm) wavelengths . Heavily hidden SF is suggested in these galaxies (e.g., 

: ITakeuchi et ^l2nnial lbl). By examining the lumi nosity functions (LFs) at UV 

and FIR from GALEX and IRAS/ Spitzei, Takeu chi. Buat. fc Burgarellal (120051 ) 
proved that the FIR LF shows much stronger evolution than that of UV, though 
^ I both evolve very strongly. This indicates that the fraction of hidden SF rapidly 

. increases toward higher redshifts up to z ~ 1. There is another important 

observable closely related to the dust emission fr om galaxies: the cosmic IR 
background (CIB). Recentlv. ITakeuchi et al.l ()2006l ) constructed the IR spectral 
energy distribution (SED) of the Local Universe. The energy emitted in the IR 
is 25-30% of the total energy budget. In contrast, the IR (from near/mid-IR 
to millimeter) contribution is ro ughly (or even more than) a half in the CIB 
spectrum (e.g.. iDole et al.ll200^ . This also suggests a strong evolution of the 
IR contribution to the cosmic SED in the Universe. Thus, understanding the 
radiative physics of dust is a fundamental task to have an unbiased view of the 
cosmic SF history. 

In general, the IR observational data are obtained from photometric mea- 
surements. Then, estimating the total IR luminosity from fluxes at discrete 
photometric bands is a n essential step to discu ss the energy budget released in 
the UV and FIR. Since FT^keuchi et al.l ( 2005al ) demonstrated good linear rela- 
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tions between MIR and total IR luminosities for several orders of magnitudes in 
luminosity. We may expect to have a good estimator, at least with a precision 
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with a factor of five. However, as well known, this is not a trivial task. Various 
estimators have been proposed based on IRAS bands, including the classical FIR 
(using 60 and 10 urn: .Helou et al . 1988), Dale et al.'s TIR (60 and 100 um: 
Dale et al. I200ltf ). and its revised version (25, 60, a nd 100 um: iDale fc Heloul 
20021'! ■ and Sa nders & Mirabel's IR (all IRAS bands: I Sanders fc Mirabe]lll99(C 
Takeuchi et al. (2005 a) examined these four estimators using a galaxy sample 
with known SEDs, and showed that they work well for normal galaxies. How- 
ever, there are a few categories of galaxies for which the FIR-based ones (the 
former two) do not give a good estimate. We focus on one of these galaxies, i.e., 
actively star-forming galaxies, and try to model their SEDs. 

Throughout this paper, we use a cosmological parameter set of {h, Qq, Aq) = 
(0.7,0.3,0.7), where h = Hq/WO [kms^^Mpc"^]. 



2. SED Model for Forming Galaxies 

2.1. Species and size distribution of dust grains produced by SNe II 

iNozawa et al.l ()2003l . hereafter N03) investigated the formation of dust grains 
in the ejecta of Population III SNe (SNe II and PISNe, whose progenitors are 
initially metal- free), taking into account the following aspects: (i) the time evo- 
lution of gas temperature is calculated by solving the radiative transfer equation 
including the energy deposition of radioactive elements, (ii) the radial density 
profile of various metals is properly considered, and (iii) unmixed and uniformly 
mixed cases in the He core are considered. In the unmixed case, the original 
onion-like structure of elements is preserved, and in the mixed case, all the el- 
ements are uniformly mixed in the helium core. [Takeuchi. Buat. & Burgarell^ 
(200^ showed that the unmixed dust production is preferred to reproduce the 



SED of a local starbursting dwarf galaxy SBS 0335—052 (discussed below). In 
addition, Hirashita et al. (,20Q5) proved that the unmixed scenario can also re- 
produce the extinction curve of a high-redshift quasar. Hence, in the following, 
we only discuss the unmixed case. The size of the grains spans a range of 
three orders of magnitude, depending on the grain species. The size spectrum 
summed up over all the grain species has a very broad distribution, and very 
roughly speaking, it might be approximated by a power law. 

2.2. Star formation, chemical evolution, and dust production 

For constructing the chemical evolution model of a young galaxy, we adopt the 
following assumptions: 

1. We use a closed-box model, i.e., we neglect an infall and outflow of gas in 
the scale of a star-forming region. 



2. For t he initial mass function (IMF), we adopt the Salpeter IMF ( Saloeteil 



^5^: (j){m) oc m with mass range of (mi,mu) = (0.1 Mq, 100 Mq). 

3. We neglect the contribution of SNe la and winds from low-mass evolved 
stars to the formation of dust, because we consider the timescale younger 
than 10^ yr. 
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4. The interstellar medium is treated as one zone, and the growth of dust 
grains by accretion is neglected. Within the short timescale considered 
here, it can be assumed safely. 

5. We also neglect the destruction of d ust grains within the young age con- 
sidered (see e.g.. lJones. Tielens. Sz H ollenba ch.1996. ). 

6. We assumed a constant SFR for simplicity. 

Using these assumptions, we calculate the chemical evolution. 
2.3. SED construction 

In this subsection, we present the constr uction of the SEP from dust. All the de- 
tails of the calculations are presented in lTakeuchi et al.1 (l2nn.'^ . lTakeuchi fc Ishiil 
(I2nn4bl 'l. and .Takeuchi et ah (2005a ). (hereafter T03, T04 and T05, respec- 
tively) . 

1. Stochastic heating of very small grains 

Very small grains cannot establish thermal equilibrium with the ambient 
radiation field. This is called stochastic heatin g. To treat the eff ect, we 
applied a multidimensional Debye model (e.g.. iDraine Sz Lill200lh to the 
specific heats of the grain species. 

2. Emission 

The em ission from dust is calculated basically according to lDraine &: AndersonI 
(|l985l ) (see T03). Total dust emission is obtained as a superposition of the 
emission from each grain species. We constructed Q{a, A) of each grain 
species from available experimental data via Mie theory. 

3. Extinction 

Self-absorption in the MIR for a very optically thick case is treated by a 
thin shell approximation (see T04). 

3. Results 

3.1. Evolution of infrared SED 

We first show the evolution of the IR SED of forming galaxies. For these 
calculation we adopted the star formation rate SFR = IMQyr^^. We adopt 
'"SF = 30 pc and 100 pc. These values are relevant when describing 'dwarf-like' 
young galaxies. The results are presented in Figure 0^. We calculated the evo- 
lution of the SED in the age range of lO'^-^-lO'^ yr. After lO'^-^^ yr, the N-MIR 
continuum is extinguished by the self-absorption in the case of rsF = 30 pc. In 
contrast, the self-absorption is not significant for rsF = 100 pc. The SEDs have 
their peaks at a wavelength A ~ 20-30 fim, which is much shorter than those of 
dusty giant galaxies at z = 1-3 detected by SCUBA. 



^The SED data are available at |http: //www.kwasan.kyoto-u. ac . jp/~takeuchi/dust .htmll 
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Figure 1. The evolution of the infrared (IR) spectral energy distribution 
(SED) of a very young galaxy. The sizes of the star forming region, rsp, is 
30 pc (solid lines) and 100 pc (dotted lines). 



4. Discussion 



4.1. Nearby forming dwarf galaxy SBS 0335—052 

It is still a difficult task to observe galaxies in their very first phase of the 
SF, especially to detect their dust emission directly. Since a recent ob serva- 
tion of SBS 0335—052 by Spitzer has been reported ( Houck et al.ll20n1 ). it is 



timely to reconsider these 'textbook objects' with new datasets. In addition, 
understanding the SEDs of these objects will shed light to the physics of inter- 
stellar matter and radiati on of high- z galaxies als o via empirical studies (e.g., 
iTakeuchi. Yoshikawa. Sz Ishii 2003 : Takeuchi et al.i.2005b ). 

SBS 0335-052 is a loca l galaxy (~ 54 Mpc) with SFR = 1.7 Moyr"! 
(iHunt. Vanzi. fc Thuanll200lh and extremely low metallicity Z = 1/41 Zq. This 
galaxy is known to have an unusual IR SED and strong flux at N-MIR. It has 
a very youn g starburst (age < 5 Myr) without significant underlying old stellar 
population (|Vanzi et al.f l200m. T03 have modeled the SED of SBS 0335-052 
and repo rted a good a greern ent with the available observations at that time. 
However, iHouck et al.l ((2004 ) presented new data of the MIR SED by Spitzer, 
and reported a deviation of the model by a factor of two or three. Their ob- 
servation indicated that SBS 0335—052 has even more FIR-deficient SED than 
ever thought. Hence, it is interesting to examine whether our present model can 
reproduce the extreme SED of this galaxy. 

We show the model SEDs for SBS 0335—052 in Figure |21 We have calcu- 
lated the SED for rsF = 10, 20, and 30 pc. The SFR is fixed to be 1.7 M© yr'^, 
and the age is 10^'^ yr. Details of the observational data are found in T03 
and T05. In the FIR regime, our model SEDs are consistent with the strong 
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Figure 2. The model for the SED of a nearby star-forming dwarf galaxy 
SBS 0335-052. 



constraint given bv lHouck et alJ ( 20041 ^ . At MIR, though we cannot give an ex- 



cellent fit to the observed data, the model SEDs roughly agree with them. The 
very strong N-MIR continuum of SBS 0335—052 is well reproduced. The dust 
mass calculated by our model at this age of SBS 033 5-052 is 1-2 x 10^ Mq, 
consistent with the observationally estimated value bv lDale et al. |(|200l3). We 
note that the mass estimation is strongly dependent on the assumed dust species 
and their emissivities, and grain size distribution. As we discussed in Section 1^ 
the continuum radiation in the N-MIR is dominated by stochastically heated 
dust emission, which is completely different from modified blackbody. There- 
fore, when we try to estimate the dust mass, we must take care to determine 
the corresponding grain properties, i.e., radiative processes and grain which are 
related to the observed SED of galaxies. 

4.2. Ly man-break galaxies 

Even in LB Gs, there is clear evidence th at they contain non-negligible amount 
of dust (e.g., lAdelberger fc Steidelll200Cir). A high dust temperature (> 70 K) is 
suggested by subsequent studies fe.g.. lOuchi et al 1 ll99fll : IChanman et al. I l200nl : 
ISawicki ^OOL l. To make a consistent picture of the dust emission from LBGs, we 
investigate the expected appearance of the LBGs with an improved dust grain 
formation of N03. 

We set the input parameters of the SED model for LBGs as follows. The 
SFR of LBGs spreads over the range of S FR ~ 1-300 M0yr~^ with a median 
of SFR ~ 2OM0yr-i (e.g., lErb et al.ll2(i03i ). Thus, the basic framework of the 
T03 model is also valid for LBGs. In this work, we consider the moderate case of 
SFR = 30 M0yr~^ over the age of lO^'^-lO^ yr. The most important information 
to calculate the IR SED is the effective size of the star forming region, but it 
is the most uncertain quantity (see T04). Since the mea n half-light radiu s of 
LBGs is estimated to be ~ 1.6 kpc from JfST observations (|Erb et al.ll2003r ). we 



use the galaxy radius as the radius of a star- forming region, and set rsF = 2 kpc 
according to T04. 
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Figure 3. The expected flux densities of an LBG. In this figure we set the 
galaxy burst age of i = 10^ yr. 



We show the observed IR/submm SEDs of LBGs at z = 2, 3, and 4 in 
Figure 01 For simphcity we only show the SED with the age of 10^ yr. The 
thick black short horizontal lines indicate the 3-0" detection limits for 8-hour 
observation by ALMA (Atacama Large Millimeter Array). Here we assumed 64 
antennas and three wavelength bands, 450, 850, and 1080 ^um. We also show 
the 3-cr source confusion limit of Herschel at 75, 160, 250, and 350 /xm bands by 
thi ck gray horizontal lines. These limits are based on 'the photometr ic criterion' 
of (iLagache. Dole, fc Pugetll2003t ) (see also iTakeuchi fc Ishiill2004al 'l . 



As discussed in T04, the detectability of LBGs is not strongly dependent on 
their redshifts. Detection at 350 /im seems impossible for moderate-SFR LBGs. 
However at longer wavelengths, if the age ^ 10® yr and SFR ^ 30MQyr~^, 
LBGs can be detected at a wide range of redshifts in the submm by ALMA 
deep survey. In the FIR, Herschel will detect the dust emission from LBGs at 
z ~ 2, but difficult at higher-z. 

4.3. Toward higher redshifts 

Based on the above discussions, we give a brief consideration on the observation 
of very high-z galaxies {z ^ 5) here. In modern hierarchical structure formation 
scenarios, it would be more reasonable to assume a small, subgalactic clump 
as a first forming galaxy. Consider a dark halo of mass ~ 10^ Mq, then it is 
expected to contain a gas with mass ~ 10® Mq. If gas collapses on the free-fall 
timescale with an efficie ncy of esp (we assume es p = 0.1), we obtain the following 
evaluation of the SFR ((Hirashita fc Huntll2004h : 



SFR ..0.1 f^) f^Y" f^y'/' [^^y^,-t] . 

VO.lAlO^Mo; VlOOpcy ^ J v; 

If Mgas ^ 10® Mq, we have SFR ~ 3(rsF/100 pc)-3/2 Mq yr"!. In addition, 
an e xtremely high-z g alaxy observed by HST has a very compact morphol- 
ogy ( Kneib et al.l 12004 ^. We also mention that, from a theoretical side, high-z 



galaxies are suggeste d to be dense and compact compared to nearby galax- 
ies (e.g., Hirashita fc Ferraral 1200^ . Thus, we consider a dwarf galaxy with 
SFR = 10 M0yr~^ as an example, and we adopt rsp = 30 pc and 100 pc. The 
age is set to be lO''^ yr. If the age is older, they will become easier to detect if 
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a constant SFR takes place. We show the expected SEDs for such galaxies at 
z = 5, 10, and 20 in Figure |1J As expected, it seems almost impossible to detect 
such objects by Herschel or ALMA. 




z=10.0 



2=20.0 





1000 10000 
glh X [^m] 



100 1000 10000 

Observed wavelength X [^^m] 



100 1000 10000 

Observed wavelength X [/^m] 



Figure 4. The expected flux densities of a dwarf star-forming galaxy located 
at z = 5, 10, and 20. Solid lines represent the SEDs for rsp — 30 pc and dotted 
lines for rsp — 100 pc. 



There is, however, a hope to observe such a small forming galactic clump 
directly: gravitational lensing works very well as a natural huge telescope. If we 
assume a lens magnification factor of 40, such a small galaxy becomes detectable 
(thin lines in Figure Since the expected SED of such a compact dwarf 
galaxy has a strong MIR continuum at their rest frame, it can be feasible to 
detect at the FIR in the observed frame. A cooled FIR space telescope is more 
suitable for such observations, and this will be a strong scientific motivation 
for a future project like SPICA. At higher-z, they can be detectable by ALMA 
survey. Even at z ~ 20, they can be detected by a standard 8-hour survey 
of ALMA, if a lensing takes place. For a practical use, we must estimate how 
frequently such lensing events occur for high-z objects. Suppose a cluster of 
galaxies at zi ~ 0.1-0.2 whose dynamical mass Mdyn is 5 x 10^^ Mq and whose 
mass distribution obeys the singular isothermal sphere. We denote the strong 
lensing cross section, i.e., the area of the region in the source plan e for which the 
resulting magnification by a cluster is larger than as a{> ji). iPerrotta et al.l 
(|2002 ) presented o"(> /x) as a function of Mdyn for ziens = 1-0. Since a{> //) oc 
{D\s is the angular-diameter distance between the lens and the source), 
we can convert their result to our condition and obtain a{> 10) ~ 30 arcsec^ 
on the source plane. This result is almost independent of the source redshifts. 
Setting the limiting flux d ensity 8^ = 1 /uJy and using the number counts of 
iHirashita fc Ferraral (120021 ') for galaxies at z > 5, we have an expected number 
of galaxies suffering a strong lensing to be ~ 1-3. Thus, we expect at least a 
few strongly lensed IR galaxies to this survey depth. 



5. Conclusion 



With the aid of a new physical model of dust production by SNe developed by 
N03, we constructed a model of dust emission from a very young galaxies. 
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The SED of a local starbursting dwarf galaxy, SBS 0335—052, was calcu- 
lated. Our present model SED naturally reproduces the strong N-MIR contin- 
uum and the lack of cold FIR emission of SBS 0335—052. Then we calculated the 
evolution of the SED of LBGs. Finally, we considered the observations of forming 
galaxies at z ^ 5. For small forming galaxies with a gas mass Mgas — 10^ Mq, 
it is almost impossible to detect their intrinsic flux by ALMA or Herschel. How- 
ever, the gravitational lensing is found to be a very effective tool to detect such 
small star-forming galaxies at z ^ 5. 
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